Comment on "Photonic Band Gaps: Noncommuting Limits and the 'Acoustic Band' 
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This Comment concerns the low-frequency dielectric 
constant e e // of 2D photonic crystals (PC's), constituted 
of circular cylinders (in vacuum) characterized by an ex- 
tremely high or infinite value of the dielectric constants. 
In Ref. and in two related papers |^,|| it was claimed 
that the quasistatic approach, based on the determina- 
tion of the slope of the " acoustic" photonic band, leads 
to a value of e e ft that is substantially smaller than the 
e e ff obtained by an electrostatic method of "homoge- 
nization" . Here we point out that this conclusion is in- 
compatible with Electrodynamics and demonstrate that 
the two approaches do give the same results, essentially 
(within computational error). 

Recently we have reported general formulas for the 
principal dielectric constants of an arbitrary 2D PC 
Q. These formulas are obtained by taking the limit 
u>, k — > directly in the wave equation for periodic me- 
dia. For isotropy in the plane of periodicity A xy = and 
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A yy , so that ei = 62 


Also, for e 


rj = 1 — / and 77(G) = —2fJ\{Ga)/Ga, where / is the 
filling fraction of the cylinders (radius a) , J\ is a Bcsscl 
function, and G is the reciprocal lattice vector. Then, 
for any isotropic 2D lattice, it follows from Ref. [[| that 
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where the exponent " — 1" implies matrix inversion. 

We have applied Eq. (1) to a square lattice; the last 
column of Table I gives the corresponding index of re- 
fraction N e ff = y/eeff- We see that there is excellent 
agreement with the static calculations ||, however poor 
accord with the results deduced from the slope of the 
acoustic band 0J^] . The authors of Ref. [JD attribute the 
discrepancy between the 2nd and 3rd columns to non- 
commuting limits e — > 00 and k — ► (k is the Bloch 
vector). However, in our approach, these limits do com- 
mute, i.e. Eq. (1) and the results in the 4th column of 
the Table I are independent on the order of taking the 
limits. Precisely the commutation of these two limits 

TABLE I. Effective indeces of refraction AT e // for three 
values of the filling fraction /. 
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provides the possibility to obtain the solution of an elec- 
trostatic problem for conductors from the solution for 
dielectics bodies of the same shape by taking the limit 
e — ► 00 [0] . We believe that the quasistatic results of Ref. 
[|lLQ are erroneous for very large and for infinite values 
of e, and that the reason for this resides in the trunca- 
tion of the generalized Rayleigh identity to finite values 
of the order I (I < 20). If this truncation is performed 
prior to taking the limit e — > 00, the stability of the nu- 
merical results does not guarantee that they are correct. 
(Indeed in Eq. ([!]) the truncation of the G values is the 
last step) .We expect that for e = 10 10 l max should 
be extremely large and that for e = 00 also Z max — > 00. 
This is based on the fact that, the higher the dielectric 
contrast between the constituents, the greater the com- 
putational effort required. On the other hand, in Ref. Q 
the contrast is given in terms of the reciprocal dielectric 
constants, leading to excellent convergence. 

The physical explanation of the disagreement Jl]-||] has 
been attributed to the extremely small wavelength within 
the cylinders. However, this argument is invalid, because 
the homogenization (in periodic systems) implies aver- 
aging over scales much larger then the Bloch wavelength 
2n/k. Local wavelength in the material is meaningless. 
The reassuring conclusion is that Electrostatics is still the 
u) — > limit of Electrodynamics - even for PC's composed 
of cylinders with an infinite dielectric constant. 
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